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Heavy-quark pair production at two loops in QCD* 
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We present updated predictions for the total cross section of top-quark pair production at Tevatron and LHC. 
For the LHC we also provide results at ^/s = 10 TeV, in view of the anticipated run in 2008 and quote numbers 
for the production of new heavy-quark pairs with mass in the range 0.5 - 2 TeV. Our two-loop results incorporate 
all logarithmically enhanced terms near threshold including Coulomb corrections as well as the exact dependence 
on the renormalization and factorization scale through next-to-next-to-leading order in QCD. 



1. Introduction 

Research on top-quark physics at hadron collid- 
ers has received great interest in the past years in 
view of the steadily improving measurements at 
Tevatron and the upcoming LHC (see Rcf. [1] for 
a recent review). In this respect, the total cross 
section for top-quark pair production is a quan- 
tity of great importance for experimental analy- 
ses and even allows for measurements of the top- 
quark mass. 

Moreover, on the theory side, the total cross 
section has been subject to numerous studies the 
motivation being improved predictions beyond 
the long-known next-to-leading order (NLO) cor- 
rections in QCD [2-4]. Recent work in this direc- 
tion has aimed at completing the next-to-next- 
to-leading order (NNLO) QCD predictions [5-9], 
at resumming large Sudakov logarithms to next- 
to-next-to-leading logarithmic accuracy [10] and, 
at estimating bound state effects [11]. Also our 
knowledge on the parton distribution functions 
(PDFs) and the precision of the top-quark mass 
determination has continuously improved over 
the last years. 

In order to study the impact of the vari- 
ous improvements on Tevatron and LHC predic- 
tions we build on the recent results of Ref. [10]. 
These approximate NNLO results for the total 
cross section are based on the complete loga- 
rithmic dependence on the heavy quark veloc- 
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ity (3 = a/I — Arn^/s near threshold s ~ Arn^. 
Moreover, they include the complete two-loop 
Coulomb corrections as well as the exact depen- 
dence on the renormalization and factorization 
scale at NNLO [12]. 

Recently, similar studies have appeared in 
Refs. [13-15]. While Ref. [13] largely follows our 
approach [10] to describe the total top-quark pair 
cross section at NNLO, Ref. [14] has limited itself 
to updating older predictions based on threshold 
resummation to next-to-leading logarithmic accu- 
racy only. Thus, Ref. [14] necessarily arrives at 
larger theoretical uncertainties. The interesting 
study of Ref. [15] on the other hand applied con- 
sistently predictions to NLO accuracy in QCD. In 
doing so, it has investigated correlations of rates 
for top-quark pair production with many other 
cross sections at LHC to quantify a potential sen- 
sitivity to the gluon luminosity. 

2. Total cross section 

The total hadronic cross section for top-quark 
pair production depends on the hadronic center- 
of-mass energy squared s and the top-quark mass 
rrit. It is given by 

cr(s, mt^) = Y^ f,/p (fif) (K) fj/p (Hf) 
i,j=q,g,9 
(g) a{mt^,fif,^i^), (1) 

where fi/p are the PDFs of the proton. The par- 
tonic cross section is given by a and (g) denotes 
the standard convolution (see e.g. Ref. [10]). 



The generally adopted procedure to estimate 
the theoretical uncertainty for a in Eq. ([T]) ex- 
ploits the residual dependence on the renormal- 
ization and factorization scale, [i^ and /Xp which 
are identified throughout this article (i.e. /x^ = 
/ijT = /i). The NLO QCD corrections for the par- 
ton cross section g and the PDFs /j/p provide 
the first instance where a meaningful error can 
be determined in this way. We define the range 
as 

cr(^ = 2TOt) - Ao-pDFifJ- = 2mf) < a{fi) 

< a(p ^ mt/2) + AapDFii^ = mt/2) , (2) 

where Acrp/^i? is computed from the variation of 
the cross section with respect to the parameters 
of the global fit (see e.g. Refs. [15-17]). 

In this contribution we employ the approximate 
NNLO result [10] to predict cross sections dH) 
and the associated uncertainty ranges ([2]) at Teva- 
tron and LHC. Let us therefore briefly comment 
on the anticipated accuracy. Our cross section 
'^NNLO (approx) ^^'^^^ ^^^^S all logarithmically 
enhanced terms In /3, fc = 1, . . . , 4 as well as the 
complete Coulomb corrections (~ 1//3. l//3^) at 
two loops for the dominant parton channels qq 
and gg and adds them on top of the exact NLO 
predictions. In this way, our predictions rely on 
exact expressions in the region of phase space 
s ~ 4?7i^, where perturbative corrections receive 
the largest weight from the convolution with the 
parton luminosities, cf. Eq. ([1]). The effect of new 
parton channels opening at NNLO {qq and qq) is 
expected to be small, cf. the qg and qg channels 
at NLO. 

The region of large energies s 3> ^m^ on 
the other hand is inaccessible within our ap- 
proach [10]. However, it is expected to give only 
small contributions in a full NNLO calculation in 
line with the observed small corrections for top- 
quark pair production together with an additional 
jet at NLO [5] which are part of the full NNLO 
correction for top-quark pair production. 

Moreover, we have also included the exact ^^ 
and /i^ scale dependence at NNLO [12] which 
can be constructed using renormalization group 
methods. For the time being, we have chosen a 
common value // for the scales, and we will ad- 
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Figure 1. The NLO and NNLO (approx) QCD 
prediction for the ti total cross section at Tevatron 
for ^/s = 1.96 TeV. The bands denote the total 
uncertainty from PDF and scale variations for the 
MRST06nnIo set [16] according to Eq. ^. 



dress the independent variation of /i^ and /if in 
a future publication. However, based on prelimi- 
nary studies we do not expect large modifications 
here. In summary, we have accounted for all nu- 
merically dominant contributions and are confi- 
dent that this provides a very good approxima- 
tion to the unknown full NNLO result as experi- 
ence from other reactions, e.g. Higgs-production 
in gluon fusion [18] shows. 

Let us next present our results for Tevatron 
and LHC. In Figs. [1] and [5] we plot the un- 
certainty range ^ comparing NLO and NNLO 
accuracy. At Tevatron (Fig. [1]) the central 

value at NNLO increases typically by 8% with 
respect to NLO. The residual scale dependence 
of '''NNLO f approx) ^^ '^'^' '^hi^h corresponds to 
a reduction by a factor of two compared to NLO. 
The overall uncertainty according to Eq. ((2]) is 
at NNLO (approx) about 8% for the CTEQ6.6 
and 6% for the MRST06nnlo PDF set. At LHC 
(Fig. [2D our ^NNLO (approx) ^^ads only to a 
small shift of a few percent in the central value 
and the NNLO (approx) band is about 6% for 
CTEQ6.6 and about 4% for MRST06nnlo, which 
exhibits again a drastic reduction of the scale un- 
certainty as compared to the prediction based on 
NLO QCD. 
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Figure 2. Same as Fig.[T]for LHC with ^s = 14 TeV. 



For phenomenological applications, the results 
of Eqs. ([H), dH) are best presented by means of 
simple formulae for the mass dependence of the 
total cross section. To that end we make the 
ansatz following Ref. [14] 

a{mt) ^ a + bx + ex + dx + ex , (3) 

where x = {rrit/GcV — 171). The parameters 
a, b, c, d, e are fitted to reproduce a in the mass 
range 150 GeV < mt < 190 GeV with a typ- 
ical accuracy of better than 0.1 per mille. For 
Tevatron and LHC the respective results for var- 
ious PDF sets are given in Tabs. [lH31 Note, that 
Eq. ^ uses a polynomial of degree four and also 
determines the parameter a for the central value 
a{mt = 171 GeV) from the fit. 

Finally, we briefly quote some NNLO (approx) 
rates for the pair-production of new heavy quarks 
in the fundamental representation of the color 
SU{3) gauge group at LHC with Vs = 14 TeV 
(see also Ref. [14]). Such particles with a mass 
rriT appear in certain extensions of the Standard 
Model and we focus on a production model which 
is entirely dominated by QCD effects. Thus, our 
cross section ct^nlo (approx) Provides a mean- 
ingful and accurate prediction because its numer- 
ical values arises largely from the threshold region 
where the logarithms In ' (3 dominate. 

In Tabs, m Owe quote the corresponding num- 
bers in the mass range 0.5 TeV < rriT < 2 TeV 
(see Ref. [14] for results to NLO accuracy). We 
observe that the scale dependence at NNLO ac- 



curacy is rather small, showing the expected good 
stability of the perturbative prediction. The rel- 
ative variation of a with respect to the PDFs, 
though, is dominating by far. Note there is 
the usual factor of two between the PDF un- 
certainty quoted by MRST06nnlo [16] and the 
CTEQ6.5 [17] PDF sets due to the definition of 
the tolerance criteria in the respective fits. The 
reason for the large observed PDF uncertainty 
is the gluon PDF being poorly constrained in 
the relevant region of large momentum fraction 
X ~ 0.1. ..0.3. This is a fact well-known to in- 
fluence many searches for high-mass particles in 
gluon fusion channels (see e.g. Ref. [15] for the 
correlation of top-quark pair production rate with 
the high mass Higgs cross section). 

3. Conclusion 

We have presented updated predictions for 
cross sections of top-quark pair production based 
on the (approximate) NNLO results of Ref. [10]. 
These represent the best present estimates for 
hadro-production of top-quark pairs, both at 
Tevatron and LHC. We have argued that the 
neglected contributions (i.e. power suppressed 
terms away from threshold and new parton chan- 
nels) are numerically small. We have found good 
convergence properties of the higher order correc- 
tions and greatly improved stability of the total 
cross section with respect to scale variations by 
our NNLO (approx) result. For applications, we 
have presented simple formulae ([3]) with 0.1 per 
mille accuracy for the mass dependence of the to- 
tal cross section in the range 150 GeV < nrit < 
190 GeV. Finally, we have applied our results to 
estimate the pair-production rates of new quarks 
heavier than the top-quark in the range up to 
2 TeV. 

The resuhs of Tabs. [TH3] for the fit of the 
mass dependence of a have also been coded in 
a C-program, which is available from the authors 
upon request. 
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Tevatron y/s = 1.96TeV 


a[pb] 


b[pb] 


c[pb] xlO^ 


d[pb] xlO^ 


e[pb] xlO^ 


CTEQ6.5 


a{n = nit) 


7.93923 


- 0.247233 


0.43143 


- 6.09398 


7.37166 




a{n = nit) - AapDF 


7.49359 


- 0.231436 


0.39979 


- 5.59529 


6.73356 




a(/i = mt) + AapDF 


8.38488 


- 0.263031 


0.46307 


- 6.59268 


8.00976 




a{n = 2mt) 


7.6572 


- 0.239371 


0.418793 


- 5.92369 


7.16791 




a{n = mt/2) 


8.00918 


- 0.247709 


0.428962 


- 6.00414 


7.2013 


CTEQ6.6 


a{^j, = mt) 


7.80984 


- 0.243547 


0.425424 


- 6.01643 


7.29114 




a{fi = nit) - AapDF 


7.36172 


- 0.22749 


0.39309 


- 5.5075 


6.6412 




a{fj, ^ mt) + AapDF 


8.25797 


- 0.259604 


0.457757 


- 6.52536 


7.94108 




a{ii = 2mt) 


7.53093 


- 0.235746 


0.412856 


- 5.84596 


7.08203 




a{fi = mt/2) 


7.87933 


- 0.244019 


0.422956 


- 5.92653 


7.11881 


MRST06nnlo 


a{^ = mt) 


8.24268 


- 0.262604 


0.46727 


- 6.70602 


8.19512 




a(/x = mt) - AapDF 


8.02331 


- 0.254847 


0.451872 


- 6.46408 


7.88561 




(j{n = mt) + AapDF 


8.46205 


- 0.27036 


0.482669 


- 6.94796 


8.50462 




a{^i — 2mt) 


7.90318 


- 0.252688 


0.450793 


- 6.47936 


7.91184 




a{n = mt/2) 


8.32601 


- 0.263279 


0.464701 


- 6.60993 


8.01315 



Table 1 

The coefficients of the parameterization ([3]) for the cross section cttsjtsjt q I'annrox") °^ ^®^' I^'^l ^^ P^ ** Tevatron 
(Vs = 1.96 TeV) using the PDF sets CTEQ6.5 [17], CTEQ6.6 [15] and MRST06nnlo set [16]. 
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LHC Vs = lOTeV 


a[pb] 


b[pb] 


c[pb] 


d[pb] xlO^ 


e[pb] xlO^ 


CTEQ6.5 


a{fj. = mt) 

a{fi ^ mt) - AapDF 

a{fi = mt) + AapDF 

cr(^ = 2mt) 

ai^i - mt/2) 


418.588 

399.508 

437.668 

413.69 

400.724 


- 11.8988 

- 11.4849 

- 12.3126 

- 11.7539 

- 11.3609 


0.19815 

0.192657 

0.203644 

0.195659 

0.188716 


- 0.267839 

- 0.261487 

- 0.274191 

- 0.264393 

- 0.254368 


3.07947 
3.01046 
3.14848 
3.03946 
2.91524 


CTEQ6.6 


(t{^ = mt) 

(j{^ = mt) - AcTpDF 

a{fi = mt) + AapDF 

a{fi — 2mt) 

a{^ = mt/2) 


419.062 

400.4 

438.121 

413.944 

401.139 


- 11.9351 

- 11.5275 

- 12.3427 

- 11.7839 

- 11.3946 


0.199035 
0.193785 
0.204284 
0.196439 
0.18954 


- 0.269327 

- 0.263554 

- 0.275101 

- 0.265736 

- 0.255779 


3.09784 
3.03731 
3.15837 
3.05688 
2.93453 


MRST06nnlo 


cr(^ = mt) 

a{^ = mt) - AapDF 

a{^ = mt) + AapDF 

(j{n = 2mt) 

a{^ = mt/2) 


449.93 

441.788 

458.073 

443.176 

431.685 


- 12.6551 

- 12.4928 

- 12.8175 

- 12.4618 

- 12.111 


0.208761 
0.206866 
0.210656 
0.205521 
0.199286 


- 0.279652 

- 0.277817 

- 0.281487 

- 0.275247 

- 0.266217 


3.1917 

3.17472 

3.20868 

3.14118 

3.03056 



Table 2 



Same as in Tab. [T] for the cross section ""mmt q faDorox") ^* LHC with start-up energy ^/s = 10 TeV. 



LHC Vs = 14TeV 


a[pb] 


b[pb] 


c[pb] 


d[pb] xlO^ 


e[pb] xlO^ 


CTEQ6.5 


(t(^ = mt) 

(j(^ = mt) - AapDF 

ain = mt) + AapDF 

a{lJL = 2mt) 

(t(^ = mt/2) 


917.522 
887.489 
947.556 
908.276 
878.259 


- 24.8588 

- 24.2054 

- 25.5122 

- 24.5978 

- 23.7272 


0.398405 
0.388932 
0.407878 
0.394104 
0.379206 


- 0.520463 

- 0.507626 

- 0.5333 

- 0.514815 

- 0.493949 


5.82374 
5.66159 
5.98589 
5.76339 
5.51362 


CTEQ6.6 


cr(^ = mt) 

cr(^ = mt) - AcrpDF 

a{^ = mt) + AcrpDF 

a{^ — 2mt) 

a{fjL = mt/2) 


920.475 
891.302 
949.648 
910.798 
881.04 


- 24.9757 

- 24.3671 

- 25.5842 

- 24.7036 

- 23.8377 


0.400681 
0.392303 
0.409059 
0.396235 
0.381364 


- 0.523783 

- 0.512828 

- 0.534738 

- 0.517914 

- 0.497136 


5.85946 
5.72046 
5.99846 
5.79166 

5.54881 


MRST06nnlo 


(t{^ = mt) 

a{iJL = mt) - AcTpDF 

a{n = mt) + AupDF 

(j{n = 2mt) 

ai^i = mt/2) 


969.128 

957.674 

980.582 

956.58 

929.84 


- 25.9797 

- 25.7648 

- 26.1945 

- 25.636 

- 24.8578 


0.412622 
0.409906 
0.415337 
0.407062 
0.393723 


- 0.534539 

- 0.530883 

- 0.538195 

- 0.527262 

- 0.508485 


5.93782 
5.88547 
5.99016 
5.86177 
5.63364 



Table 3 



Same as in Tab. [T] for the cross section o"i>^nLO fanorox) ^* LHC with ^/s = 14 TeV. 





only scale uncertainty 


only pdf uncertainty 


total uncertainty 


TOT[TeV] 


min 


max 


5%] 


min 


max 


5% 


min 


max 


5[%] 


0.5 


4345. 


4472. 




4287. 


4656. 


4 


4160. 


4656. 


6 


0.6 


1561. 


1601. 




1526. 


1676. 


5 


1486. 


1676. 


6 


0.7 


634.1 


649.2 




616.0 


682.5 


5 


600.8 


682.5 


6 


0.8 


282.3 


288.5 




272.6 


304.4 


6 


266.4 


304.4 


7 


0.9 


134.5 


137.2 




129.3 


145.1 


6 


126.6 


145.1 


7 


1.0 


67.64 


68.94 




64.81 


73.08 


6 


63.50 


73.08 


7 


1.1 


35.45 


36.17 




33.93 


38.41 


6 


33.22 


38.41 


7 


1.2 


19.23 


19.65 




18.38 


20.91 


6 


17.97 


20.91 


8 


1.3 


10.74 


10.99 




10.26 


11.72 


7 


10.01 


11.72 


8 


1.4 


6.147 


6.301 




5.862 


6.741 


7 


5.708 


6.741 


8 


1.5 


3.589 


3.687 




3.417 


3.957 


7 


3.319 


3.957 


9 


1.6 


2.130 


2.192 




2.021 


2.363 


8 


1.959 


2.363 


9 


1.7 


1.282 


1.322 


2 


1.212 


1.432 


8 


1.172 


1.432 


10 


1.8 


0.781 


0.806 


2 


0.735 


0.878 


9 


0.710 


0.878 


11 


1.9 


0.480 


0.497 


2 


0.450 


0.544 


9 


0.433 


0.544 


11 


2.0 


0.298 


0.309 


2 


0.277 


0.340 


10 


0.266 


0.340 


12 



Table 4 

The NNLO (approx) cross section of Ref. [10] in fb for the pair-production of a (new) heavy quark with mass 
rriT at LHC {^/s = 14 TeV) using the MRST06nnIo PDF set [16]. 5 is the relative uncertainty with respect to 
the central value: 5 — 100 x (max — min)/(max + min). 





only scale uncertainty 


only pdf uncertainty 


total 


uncertainty 


TOT [TeV] 


min 


max 


S[%] 


min 


max 


5[%] 


min 


max 


S[%] 


0.5 


3921. 


4037. 




3639. 


4436. 


10 


3522. 


4436. 


11 


0.6 


1402. 


1440. 




1275. 


1604. 


11 


1238. 


1604. 


13 


0.7 


568.0 


582.5 




508.6 


656.5 


13 


494.1 


656.5 


14 


0.8 


252.4 


258.2 




222.5 


294.0 


14 


216.6 


294.0 


15 


0.9 


120.3 


122.8 




104.5 


141.0 


15 


102.0 


141.0 


16 


1.0 


60.48 


61.66 




51.94 


71.37 


16 


50.76 


71.37 


17 


1.1 


31.78 


32.30 




26.95 


37.64 


17 


26.44 


37.64 


17 


1.2 


17.25 


17.57 




14.50 


20.57 


17 


14.21 


20.57 


18 


1.3 


9.626 


9.802 




8.023 


11.58 


18 


7.847 


11.58 


19 


1.4 


5.503 


5.614 




4.555 


6.673 


19 


4.444 


6.673 


20 


1.5 


3.208 


3.277 




2.630 


3.925 


20 


2.560 


3.925 


21 


1.6 


1.902 


1.946 




1.545 


2.347 


21 


1.502 


2.347 


22 


1.7 


1.144 


1.173 




0.921 


1.426 


22 


0.892 


1.426 


23 


1.8 


0.696 


0.715 




0.554 


0.876 


23 


0.535 


0.876 


24 


1.9 


0.428 


0.441 




0.337 


0.544 


24 


0.324 


0.544 


25 


2.0 


0.265 


0.274 


2 


0.206 


0.342 


25 


0.198 


0.342 


27 



Table 5 
Same as in Tab.g]using the CTEQ6.5 [17] PDF set. 



